Evaporation of a black hole off of a tense brane 
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We calculate the gray-body factors for scalar, vector and graviton fields in the background of 
an exact black hole localized on a tensional 3-brane in a world with two large extra dimensions. 
Finite brane tension modifies the standard results for the case with of a black hole on a brane with 
negligible tension. For a black hole of a fixed mass, the power carried away into the bulk diminishes 
as the tension increases, because the effective Planck constant, and therefore entropy of a fixed mass 
black hole, increase. In this limit, the semiclassical description of black hole decay becomes more 
reliable. 
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A distinct signature of low scale quantum gravity mod- 
els [l| is the possibility of black hole production in high 
energy collisions Such processes could be probed 

in particle accelerator experiments in the very near fu- 
ture. Recently a great deal of effort has been expended 
for the precise determination of observational signatures 
of such events. Chief among them are the black hole 
production cross section and the Hawking radiation spec- 
trum. To date, virtually all the work has been done for 
the idealized case where the brane tension is completely 
negligible. This is because it is generically very hard to 
obtain exact solutions of higher-dimensional Einstein's 
equations describing black holes on branes with tension. 
On the other hand, one generically expects the brane ten- 
sion to be of the order of the fundamental energy scale 
in theory, being determined by the vacuum energy con- 
tributions of brane-localized matter fields. 

Recently, a metric describing a black hole located on 
a 3-brane with finite tension, embedded in locally flat 
6-dimensional (6D) spacetime was constructed In 
spherically symmetric QD Schwarzschild gauge, this met- 
ric can be written as 



dt 



sm 



sin^ (/)(rfx^ + B'^ sin^ xdip^)] } 



(1) 



The parameter B — 1 — 2ttM* measures the deficit angle 
about the axis parallel with the 3-brane, in the angu- 
lar direction ip, because the canonically normalized angle 
i/y = ip/B runs over the interval [0, 2tt/B]. Here A is the 
brane tension, and is the fundamental mass scale of 
6D gravity, defined as the coefficient of the Ricci scalar 
in the 6D action J d'^Xy^M^Re/2 e Sqd- The black 
hole horizon is at 



rh 



(2) 
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given in terms of the black hole ADM mass 
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4^ 



1 



M 



BH 



(3) 



Thus the main effect of the brane tension is to change 
the relation between the black hole mass and horizon 
radius. Then in terms of the geometric quantities, the 
metric appears the same as the QD Schwarzschild solu- 
tion. However, because of the deficit angle, at large dis- 
tances (r 3> Th) the geometry asymptotes to a conical 
bulk space. 

Following 0, Q I we take as the black hole production 
cross section its horizon area, 

a^^rl. (4) 

It is interesting to explore how this quantity depends on 
the brane tension A. To do this, one must be careful be- 
cause one must first appropriately compactify the asymp- 
totic geometry of the solution ([T]), in order to be able to 
think of it as a small black hole on a brane in large com- 
pact dimensions. In this case, the volume of the compact 
dimensions will also depend on the deficit angle B. The 
precise details depend on the compactification. 

One simple case would be to imagine ending the space 
far from the black hole ([T]) on a cylindrically symmetric 
4-brane [9| and imposing relfection symmetry about it, 
in which case the enclosed 2D volume would be V'^'^^ ~ 
2'kL'^ /B, when the black hole is much smaller than the 
radial size of the extra dimensions, rh ^ L. Because 
the 4D Planck mass M4 is determined by the Gauss law 
equation 



u 



(5) 



where rc is the conventional 61? Schwarzschild radius. 



it is then possible to vary the tension on our 3-brane, 
and, by matching that with a change of the tension on 
the cylindrical 4-brane, simultaneously hold and A/* 
fixed. This would decrease B and enhance the cross sec- 
tion for a fixed value of and Mbh , while the decrease 
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of from the reduction of the opening angle of the 
cone can be compensated by the increase in the hnear 
size of the extra dimension L, at least in principle. 

Obviously, there are restrictions to how far one might 
take these limits in order to not exceed the sub- millimeter 
bounds on the corrections to gravity, which were recently 
reviewed in [lol |. In general, the situation is even more 
complicated because one must carefully consider the de- 
tails of compactification, and its topological and geomet- 
ric aspects, as well as the presence and deployment of 
any additional branes in the bulk, to determine various 
possible contributions to the Gauss law ([5|) to see how 
they compare to the tension-dependent contribution of 
the black hole-bearing brane. These detailed considera- 
tions may also alter the asymptotic form of the metric ((T]) 
far from the black hole, as they generally are expected 
to do so in any extra-dimensional framework. This dis- 
cussion shows that the precise dependence of the black 
hole production cross-section on the brane tension is in 
general not simple, but is sensitive to the details of com- 
pactification. For the purposes of our discussion, we can 
however ignore this issue, and imagine that we can ar- 
range for a compactification where M*, M4 and A are 
freely tunable parameters. With these assumptions we 
can treat the solution ([1]) as a good approximation for a 
small black hole on a brane. 

The smoking gun of black hole production in collid- 
ers, and hence of low scale quantum gravity, would be 
their decay via Hawking (or rather Hawking-like ra- 
diation. It is therefore essential to accurately calculate 
emission spectra for fields of various spin emitted by a 
black hole, whose precise nature is encoded in the black 
hole gray-body factors. Now, from the metric ([T]) it is 
clear that the gray-body factors for the brane-localized 
fields remain formally identical to those calculated for the 
tensionless brane, with the only difference being in the 
relationship between the black hole mass and the hori- 
zon size, and therefore the geometric cross section, as 
encoded in the Eqs. ^ and (g]). 

However, the fields propagating in the bulk will have 
diff'erent gray-body factors. While in the simplest models 
with large dimensions only the fields from the graviton 
multiplet are taken to propagate in the bulk, in order to 
suppress a rapid proton decay, one needs to physically 
separate quarks from leptons [5[. The simplest way to 
realize this is to localize quarks and leptons on different 
branes in the bulk. Then, since gauge fields must interact 
with both quarks and leptons, they must also be propa- 
gate through the section of the bulk between the quark 
and lepton branes. When a small black hole resides on 
one of those branes, having been created in a collision of 
matter particles localized on them, it will emit the gauge 
fields and its possible superpartners as bulk fields. Thus 
it is important to to include bulk fields in the detailed 
black hole evolution in order to develop the tools for the 
study of more realistic braneworlds. 



Let us first consider bulk scalars. For small hot black 
holes we can ignore the scalar field masses and simply 
analyze the equation of motion of a massless bulk scalar 
□5' = 0, where 5* propagates in the background llj. We 
can separate this equation into the radial and angular 
parts. 



1 d 
dr 



[Fr'±R(r))+i^~^)Rir)^0, (6) 



and 



A4r(0,0,X» = -rym0,X,^) 



(7) 



where A4 is the Laplacian on the deformed 4-sphere 
dsl = de^ + sin^ + sin^ (l>{dx^ + sin^ x^^)] . 

Here rj is the separation constant, F = (1 — (^)^), 
and B measures the deficit angle. For tensionless brane, 
_B = 1, the angular equation reduces to the one for the 
spherically symmetric scalar field. The solutions can be 
written as the expansion in hyper-spherical harmonics 
with eigenvalues rj = L{L + 3) For each fixed L 

there is (2L + 3)(L + 2){L + l)/6 states. When B ^ 1, 
spherical symmetry is broken, and eigenvalues and degen- 
eracies will change, although the total number of states 
will remain the same. 

To proceed with the general case B 7^ 1, we sep- 
arate the Eq. ^ further, by writing Y{6,4),Xt'>P) = 
0(0)<I'((?!))r(x)S(7/'). This yields four angular equations 



1 
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^edo 



sm 
1 d 
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1 d^ _ 



r - 



-»?2r, 

-771 s. 



(8) 



The solutions will therefore be classified by four quantum 
numbers, {L,li,l2,m) which determine the four eigen- 
values (77,773,772,771) respectively. L, /i, I2 and m are 
non-negative integers (negative m's will have the same 
eigenvalues as positive ones, so we can restrict to ttj > 
without any loss of generality) , and we have the hierarchy 
L > li > I2 > \m\. In the limit oi B = 1, spherical sym- 
metry implies 77 = L{L+3), 773 = /i(/i + 2), 772 = ^2(^2 + !), 
771 = m^. In this limit spherical symmetry implies that 
the eigenvalue 77 depends only on the quantum number 
L. However when _B 7^ 1, spherical symmetry is broken 
and 77 depends on two quantum numbers, i.e. m and L. 
To determine the black hole power, we thus need to de- 
termine degeneracies iV/, „i in the states labelled by m 
and L. We find 



{L + 2){L + 1) 



(9) 



m 7^ 



Nl%^{L-m + 2){L 
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We emphasize here that only the eigenvalue rj appears 
in the radial equation However, because the eigen- 
values are coupled, to find 77 one still needs to solve the 
full system of equations ([5]). To determine 77, we adopt 
the numerical technique following the Chapter 17.4 of 
and rewrite the first three equations in ([5]) as: 



which are respectively 



VA!^'^. (13) 



1 



sm a 



da 



Here, r* is a "tortoise" coordinate defined by dr* — 
h r]k+i)P = . (10) dr/F. Choosing the boundary condition A^^^ = which 



2 

sm a 



The parameter k takes the values k — 3, 2, 1 and the 
angle a stands for ^, </> and x for the first, second and 
third equation in ([5]) respectively. After setting x = cos a 



and changing the function P to P 

n = (-(fc - 1) + v/(fc 



(1 - a;2)"T with 



1 



4?7fc)/4, we obtain 



(l-a;2)— r-(4n+fc+l)x— r+(r/-2n(2n+fc))r = , 



dx^ ^ ' ' ' dx 
with the boundary conditions 



(11) 



1 



x = 



dT _ {ri' - 2n{2n + k)) 
dx ~ (4n + A; + 1) 
dT 

T = 0, or — = . 

dx 



(12) 



The two different boundary conditions at x = corre- 
spond to symmetric (i.e. dT/dx = 0) and antisymmetric 
(i.e. T = 0) solutions. We can then use the shooting 
method to solve for the eigenvalue r;. In Fig. [1] one can 
see how this eigenvalue changes with B in nine branches 
of L. The m = branch does not change with B. How- 
ever 771 ^ branches increase quickly as the brane tension 
A increases, and so B decreases. We note here that the 
increase in the eigenvalue will result in the reduction of 
the power of emitted radiation. 




FIG. 1: The figure shows how the eigenvalue 77 changes with 
B. L changes from to 9. 



ensures that near the horizon the solution is purely in- 
going, we can numerically integrate Eq. ([6]) using the 
forth-order Runge-Kutta method. From this solution we 
can calculate the absorption ratio, 



1^ 



L.m I 



= 1 - 



A 



(00)' 



A, 



(00) 



(14) 



and then, using the principle of detailed balance, write 
down the spectrum of Hawking radiation. 



dtdoj 



27r 



(15) 



where T;, = 



is the Hawking temperature. In Fig. 



we display the results of numerical integration, which 
show the variation of the Hawking radiation spectrum 
with brane tension, while holding r/j and M^, fixed. For 
a tensionless brane, B = 1, the result coincides with the 
existing literature (c.f. in [l^l). Increasing the brane ten- 
sion, and so decreasing S, while assuming that M4 and 
AU are held fixed, reduces the emitted power. We note 
here that most of the 7n 7^ modes have large values for rj 
and are thus suppressed with respect to the tti = mode, 
in agreement with the statement that most of Hawking 
radiation from a non-rotating black hole is emitted in the 
s-wave channel [1] (see also 14 1). 




Having thus determined the angular eigenvalues 77, we 
can use Eq. ^ to find the asymptotic form of the radial 
wavefunctions in the far-field zone and near the horizon. 



FIG. 2: Hawking spectrum of a QD bulk scalar as a function 
of brane tension. B changes from B — 1 (tensionless brane) 
to B = (deficit angle equal to 2-k). 
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We should stress that our method of calculating the 
Hawking radiance by bulk mode expansion is completely 
equivalent to first performing the Kaluza-Klein reduction 
of bulk scalar perturbations and then expanding in the 
50(3) representations reflecting the black hole horizon 
symmetry. This way of thinking about organizing the 
calculation is particularly useful in computing the higher 
spin emissions, which we turn to next. 

The calculation of Hawking radiance for bulk vectors 
and tensors goes similarly to the one for scalars, with 
the exception that one has to properly account for the 
helicity multiplets. Namely, in the case of massless AD 
vectors and tensors, gauge symmetries ensure that the 
only propagating modes are, in the traceless-transverse 
Lorentz and de Bonder gauge, respectively, the helicity- 
1 and helicity-2 modes, like the usual Maxwell and Ein- 
stein fields. However, in the higher-dimensional case they 
are accompanied by the Kaluza-Klein towers of massive 
modes. These have less gauge symmetry and hence con- 
tain more propagating degrees of freedom, including ex- 
tra scalars, and scalars and vectors in the bulk vector 
and bulk tensor cases, respectively, which span the full 
massive AD vector and tensor multiplets [15i] . 

Now, in computing the Hawking radiance from small 
hot black holes we can neglect the AD mass for most of 
these modes since they lie well below the cutoff while the 
black hole's temperature is close to it, but these extra 
modes still carry away the black hole energy. Thus we 
must take them into account. A simple way to do this is 
to dimensionally decompose the relevant representations 
into AD tensors, vectors and scalars, as has been done in 
[3 [ill J with respect to the transverse space to the 
direction of motion. In practice, this means counting the 
tensor representations over the deformed transverse S"* in 
Eq. ([1]) , because by the axial symmetries of the problem 
the energy will flow radially outward from the hole. Thus 
the radial direction and the time can be factored out. 
Then, by going to the unitary gauge, we will find that 
each degree of freedom obeys a field equation similar to 
the one for scalars. Working out the details, one can show 
that the master equation for the radial part of a higher 
dimensional gauge boson or graviton field (see 
18|) propagating in the background of the black hole ([T]) 
is 



+ (u^ - Vir))^ = 0, 



where 



V{r) 



, n 2 



4(1 



(16) 



(17) 



The constant p is given by for scalars and tensor gravi- 
tons, 2 for gravi- vectors, 1/2 for gauge vectors and 3/2 for 
scalar reductions of bulk vectors. For example, as a quick 
check we can see that setting p = and = r^R in ITBl) 
we readily recover Eq. One further finds that for 



the higher spin fields the angular quantum number L is 
restricted by the spin, so that its range is 



m 17, 18, 1 



0, 1... for scalar perturbation , 
L = ^ 1,2... for vector perturbation, 
2,3... for tensor perturbation . 

Finally for gravitational scalar perturbations one finds 

Mm 



V = F 



Q{r) 



(18) 



with 



Qir) 



6A00r-^ + 1920(77 - 4)r"^ - 1920(7? ~ 4)r"^ 
+ 16(r7-4)3-|-96(r/-4), 
H{r) = (r;-4) + 10r"^ (19) 



What then remains is to determine the eigenvalues 77 
and the and the number of states for each fixed 77. This 
number depends on the helicity, and for the scalars the 
degeneracy Nf^f\ is the same as in Eq. (ITOl) . For the 
vectors, one needs to account properly for gauge fixing 
when counting the solutions of the angular equation of 
motion [l3|. In the end the vector degeneracies depend 
only on the helicity and not on how it is lifted into the 
bulk, yielding 



L(3L + 5) 



(20) 



771 7^ 



TVn, - 3(L - 777 + 2)(L - 777 + 1) . 



In the case of transverse-traceless tensors, the problem 
of solving the gauge fixed angular equation is more dif- 
ficult. To determine the degeneracies, we need to solve 
V^VaTpjy — satisfying gauge conditions V^T^j/ = 
and T'^p = 0, where V"Vq is the covariant Laplacian 
on S"*, and all the indices a,^,v run over the deformed 
transverse sphere of Eq. ([1]), as we noted above. Now, 
below some high (angular!) momentum cutoff, the to- 
tal number of states is the same as in the spherically 
symmetric limit, computed in [20| and summarized in 
the Table I of that paper. The only effect of broken 
spherical symmetry is splitting between some of the lev- 
els that were degenerate on S^. Then for i > |777| > 2 
the transverse-traceless tensor can be viewed as an array 
of five scalars, and we can simply count up the number 
of states for each of them and add them up. This yields 
Nf^ = 5(L - 771 -h2)(L- 777+1) for L > |777| > 2. What 
remains then is to account for the states with ttt, — 0, 1. 
We can do this by taking the total number of states from 
[20| . subtracting all the states with |77i| > 2, whose num- 
ber we just worked out, and recalling that the number of 
states with L = 2, 777, = is really the number of indepen- 
dent s-wave modes in the full 6D graviton multiplet at 
some fixed mass level. Since the number of independent 
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propagating modes of a 6D graviton is nine, we must fix 
N{L = 2, m = 0) = 9. This then uniquely determines 
the degeneracies to be 



m = 

m — I 
m > 2 



^tcn 

L,rn 

^tcn 

^tcn 
L.rn 



_ (L- l)(5£ + 8) 
2 

= (i- l)(5L + 6), 

= 5(i - TO + 2)(i- 771+1), 



(21) 



Then using the same asymptotic form of the radial 
wavefunctions as for the scalars (fT3|) and defining the 
absorption ratio for each helicity as I^^L^mP = 1 ~ 

^Ag^^'/A-^^^ , where j labels the hehcity, we can write 
the Hawking radiance spectrum for each mode as 



dtdw 



(22) 



For vectors, j takes values for scalar and Lorentz vec- 
tor, whereas for tensors it runs over all scalars, vectors 
and transverse-traceless tensor modes with fixed quan- 
tum numbers L and to. This is the total power emitted 
in the frequency band dw. For the power emitted per 
degree of freedom, we ought to divide the result ([22]) by 
four for a gauge boson and by nine for graviton, since 6Z3 
bulk vectors and tensors have 4 and 9 degrees of freedom, 
respectively, as we have discussed above. 

We have computed the power ([55]) numerically for the 
relevant vector and tensor modes. In Fig. [3] we show 
the dependence of the Hawking radiation spectrum for 
bulk vectors on the brane tension, while in Fig. 3] we 
display the dependence of the radiance on tension for 
the graviton field. Finally, for completeness, in Fig. [5] 
we show the Hawking spectrum of brane-localized fields, 
which for a fixed horizon size rh does not depend on the 
tension. Of course, the horizon size rh depends on the 
tension as given in Eq. ([2]) . In order to get the power per 
channel from the total power displayed in Fig. ([5]), we 
need to halve the fermion and gauge boson distribution. 

In sum, collider searches for black holes may be an ex- 
citing and interesting source of information about extra 
dimensions, complementing possible astrophysical tests 
(for recent exploration of those, see [2lj). In this note 
we have explored the Hawking decay channels for a re- 
cently constructed exact black hole localized on a 3-brane 
0. We see that brane tension may alter significantly 
the power output of small black holes located on the 
brane. While for the non-rotating black holes the domi- 
nant channels are still the brane-localized modes, in re- 
alistic models with split fermions as in the gauge 
boson contribution to black hole radiance may be signif- 
icant. Since they need to be in the bulk to communicate 
between different (families of) quarks and leptons, the 
gauge boson emissions may be a sensitive probe of deficit 
angle and therefore brane tension. Specifically, we see 
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FIG. 3: Hawking radiation spectrum for a bulk gauge boson 
as a function of the brane tension. B changes from B = 1 
(tensionless brane) to B = (deficit angle equal to 2ti). 
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FIG. 4: Hawking radiation spectrum for a bulk tensor as a 
function of the brane tension. Again, B changes from B = 1 
(tensionless brane) to i? = (deficit angle equal to 2tt). 



that the power emitted in the bulk diminishes as the ten- 
sion increases. This may appear slightly odd at first sight 
since increasing the tension for fixed M4 strengthens ef- 
fective bulk gravity However, when the horizon rh 
is held fixed, this comes about because the gravitational 
potential barrier depends on the eigenvalue 7y, which as 
seen in Fig. ([T]) increases as the parameter B decreases, 
and tension increases. Therefore the potential barrier 
grows reducing the rate of energy loss. Similarly, at fixed 
black hole mass increasing the tension would simultane- 
ously increase the horizon size, lowering the temperature 
and increasing the black hole entropy. Hence in this case 
a black hole will also radiate away its mass more slowly. 
If we ever produce small black holes at the LHC, the 
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FIG. 5: Hawking radiation of brane localized fields: a scalar, 
a fermion and a vector gauge boson. For fixed horizon size 
they are not affected by the finite tension of the brane. 



total Hawking radiance will be a sum of the power emit- 
ting along the brane and off of it. In simple toy models 
the emission into the bulk would be subleading because 
of the s-wave dominance, leading to essentially negligi- 
ble bulk effects. However in more realistic models with 
split fermions and gauge fields in the bulk, the balance 
between contributions from the brane and bulk modes 
may end up being significantly tilted in favor of the in- 
creasing number of bulk modes. While the precise details 
of course would be very model dependent, in such cases 
the off the wall contributions may be significant, and so 
the effects we have described here may play an impor- 
tant role, warranting a careful and detailed examination 
of Hawking radiation losses. 
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